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Abstract 

I discuss the universal aspects of scaling in period-doubling sequences in 
families of maps of the real line possessing non-integer degree. I show that 
the scaling behaviour in both the orbital and parameter spaces is governed 
by the same sequence of eigenvalues of the linearized renormalization oper- 
ator. These eigenvalues are smooth functions of the degree of the maximum 
of the map. 

1 Introduction 

In his seminal paper [1], Feigenbaum showed that period-doubling sequences in 
families of maps of the real line of the form 

f,{x) = X-\xf {d = 2,4) (1) 

have the asymptotic behaviour 



Afc+i -Xk = AXk ~ j^ (2) 

as A; — > oo. Here X^ is the smallest parameter value at which fx possesses a 
superstable 2^ cycle, superscript [n] indicates n-fold composition, and Gq and bo 
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are constants, (oq and 60 are of course dependent on 0?, but as I am not interested 

in their value, I suppress the dependence from the notation), a and h were claimed 

to be universal^ that is, dependent on d only. Additionally, Feigenbaum gave an 

argument ([1], section 5) that the rate of convergence of a^ = (j)k/4>k+i to its 

limit is also 5. This amounts to claiming that equation (3) is the first term in an 

expansion 

bo bi 

"^^ ~ a{dY ^ {a{d)8{d)Y ^ ■ ■ ■ ■ ^^' 

However, Feigenbaum had evidence was that this is true for d = 2 only. The 
present study arose out of a desire to understand this discrepancy, and more 
generally, the dependence on d of the spectrum of scaling exponents. The case 
d = 2 has been studied by Mao and Hu [9], Liu and Young [8] and by Reick 
[11]. Here I present results of a numerical study of the problem, and a plausibility 
argument to justify the assumptions used the numerical study. 

2 Numerical results 

We postulate the following forms for the behaviour of AX^ and (pk'. 



AA. = Et^ (5) 

8 = "«■ 

^^ = t^- (6) 

8 = "8 

To ascertain the validity of these expansions, I accurately computed from four 
to eight of the quantities 5i and a^ for about 50 values of d between 1 and 10. 
Here 5i and a^ are the constants to be determined, with Sq = 5 and ao = «• The 
exponents are named in order of increasing magnitude. I found it necessary to 
compute typically 15 to 20 superstable parameter values A/;, for which I used a 
Newton iteration method [3]. The use of high precision arithmetic (50 to 100 deci- 
mal places) becomes essential, since up to 2^° iterations of fx are being computed, 
and there is otherwise an unacceptable accumulation of roundoff errors. 

I then estimated the exponents by the method described by Mao and Hu 
[9], again using high precision techniques. I illustrate this by the case of the 5 
scaling exponents, the case of a being exactly analogous. The problem reduces to 
solving a linear system and a polynomial equation, as follows. Assuming that N 
scaling exponents are desired, 2A^ — 1 values of AA are required. Let us call these 
AAj_|_i, . . . , AAj_|_2jv. Fliminating the constants a^ from the definition of the scaling 



exponents results in 



-AA,+2 +AA,+3 
-AAj'+jv +AAi+jv+i 



■1)^AA,+Ar 

■n^AA 



i+N+1 



D^AA 



i+2N 



t2 



-AA, 
-AA,+i 



-AA,: 



i+N-l 



where ti = J2i=i ^i^ ^2 = 1^8=^? ^i^j stc. This equation is solved for {ti^t^ 
and the roots of the polynomial 



N _ ™JV-1 



h + --- + {-iy'tN 



(7) 



>^jv}, 



(S) 



are found. The N roots are the desired scaling exponents 5i. This calculation can 
also be viewed as the construction of a Pade approximant to the analytic function 
/(^) = I^i^i AAjX', and the subsequent determination of the poles 5i of f[x). This 
point of view suggests that in the case of integer 0?, where 5 is known to very high 
accuracy [4], elimination of the pole 5 by multiplication of f{x) by 1 — x/5 might 
achieve greater accuracy in the calculation of the higher poles 5i. However, I have 
not found a sufficient increase in precision to make this procedure worthwhile. 

A selection of numerical values are shown in Figure 1, and the full results are 
shown in Figures 2-5. For reasons that will become apparent, it is most useful to 
consider 5i/5 and ai/a^i = 1,2, . . ., rather than the exponents themselves. Only 
as many decimal places as are believed correct are quoted. 
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Figure 1: The dominant scaling exponents for even integer d. Here 'origin' refers 
to the explanation of the numerical values in terms of the eigenvalues /ij-, and the 
gaps are unobserved scahng exponents. 



3 Theory 

I believe that all the above results can be understood by an extension of an argu- 
ment of Feigenbaum [2]. I define the usual doubling operator T, by 

{Tf){x) = afof{a-'x), (9) 

where a is the (assumed known) value appropriate to the degree d of the maximum 
of /. A Taylor expansion in A of fx about Aqo gives 

h = /a. + (A-Aoo)|^ + --- (10) 

= go{x) + {X - Xo.)ho{x) + ■ ■ ■ . (11) 

gk and h^ may now be defined by 

{T''h){x) = gk{x) + (A - X^)hk{x) + • • • , (12) 

and it is known from the work of Feigenbaum [2] that as A; — > oo 

9k{x) -^ g{x) (13) 

hk{x) -^ {DT,){hk-i){x). (14) 

Here g{x) satisfies Tg = (/, DTg is the linearization of T about (/, and q~^ = (/(O). 
It is known that Feigenbaum's 8 is the dominant eigenvalue of T^, and governs 
the approach of Xk to Aqo- I will now compute the effect of the subdominant 
eigenvalues. 

Let us expand ho[x) in eigenvectors ej(a;), (j = 1,2,...) of DTg. Letting 
DTcj = fijCj with 6^(0) = 1 and fii = 5 gives 

oo 

ho{x) = ^ajej(x) (15) 



j=i 



a,; 



f^] \k. 



hk{x) = a,5'[e,{x) + J2^{^)'e,{x)]. (16) 

Thus as A; — > oo, 

a'flia-'x) ~ g{x) + (A, - X^)a,S'[e,{x) + ^ ^(^)'=e^.(x)], (17) 

and so 






a,; 



1^] \k-\-l 



X^-X,r.a-'S-'^[l + Y.^i^)'^]-\ (18) 

Inverting the sum in square brackets shows that 



oo 



Aoo-Afc~^6,(^^^ (19) 



Figure 2: The scaling exponents ajU s.d, z = 0, 1, 2, 3. 
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Figure 3: The scahng exponents Sivs.d^ i = 0, 1, 2, 3, 4. 



Figure 4: The ratios aj/a vs. d^i = 0, 1, 2, 3. 
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Figure 5: The ratios 5i/5 vs. o?, z = 0, 1, 2, 3, 4. 



where each 8^ is of the form 8 YYk=i f^k"^'' ■> ('^fc non-negative integers), and {6^} are 
constants which are comphcated functions of {oj}. 

This completes the computation of 8^ in terms of {/ij}. It is clear that the 
behaviour can be pictured as an approach of /aj. to /a^, which for an arbitrarily 
chosen parameterisation of /a, will not be along the eigendirection corresponding 
to 8. Nevertheless, the more general approach is still governed by the largest 
eigenvalues of DTg . 

I now perform a similar calculation for a^. Applying the doubling operator 
2^~^ times to f\^ gives 



a^ft\a-^^-'^x) ~ g{x) + (A^ - \^)a,8^-\e,{x) + ^ -±{^^f-\^{x)]. (20) 



a,- 



A'j \fcl-l 



a'ci^, ~ 1 + (A, - \^)a-'8\l + Y: ^{^f]-' (21) 



k-l 



(recalling that (f)k = /a (0))- Now performing an inversion as above gives 

oo 

where Cj are constants and 8^ are the same scaling exponents as above. The 
important result is that the one set of eigenvalues /ij governs both the parameter 
and the orbital scaling. The latter simply have an extra factor of a. 

There is a curious phenomenon in the case d = 2. The largest eigenvalue 
less than one takes the value exactly a~^ . In other words, there is an eigenvalue 
crossing at o? = 2. As expected, the effect of the eigenvalue is observed in the 
orbital scaling, but surprisingly, it is not observed in the parameter scaling. I 
conjecture that this is because the corresponding amplitude a^ is zero (or very 
small). This effect nevertheless deserves more detailed study. This appears to be 
the only case where /ij = a~^ for any j and k, though there are close approaches at 
0? = 4, 6 and 8. These close exponents make the numerical calculation difficult, and 
this explains the oscillations seen in the graphs at larger d. I believe that a more 
accurate calculation would show that all the exponents are monotonic functions 
of d. 

The eigenvalues /ij can be computed independently for integer d via a numer- 
ical approximation to DTg [4]. This can be achieved this by first computing a 
finite-dimensional approximation to g{x) in the form g{x) = J2q gix'^\ This is 
found by solving Tg — (/ = by a Newton iteration in the space of Taylor coeffi- 
cients gi. Similarly DTg is represented by a A^ X A^ matrix, and all its eigenvalues 
computed. The fact that an infinite-dimensional problem has been truncated to 
finite dimension means that some extra eigenvalues are introduced, but the rele- 
vant eigenvalues are usually readily identified by comparison with the bifurcation 
data. It is known that a~^ is an eigenvalue for integer k [1], and Feigenbaum 
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conjectured that these are all the eigenvalues of modulus < 1. It is clear now that 
this is not the case. The eigenfunctions corresponding to even integer k are even, 
so that for odd integer o?, our expansions cannot represent the eigenfunctions, and 
the eigenvalues are not found. However, this is not a difficulty as these eigenvalues 
are not relevant to the bifurcation data. 

I list a few of the largest eigenvalues fii for even integer d at the bottom of the 
tables in Figure 1. Observe that eigenvalues not of the form a~^ are also present. 
This then is the explanation of the the discrepancy between the behaviour at o? = 2 
and 0? = 4 mentioned in the introduction. 

Of interest is the question of the limits o? — > 1 and o? — > oo of the scaling 
spectrum. Collet et al [6] have shown that (^ — > 2 as o? — > 1. Eckmann and 
Wittwer [7] have shown that 8 — >sa 30 as (^ — > oo. On the question of this limit, 
see also [12]. It is also known [10] that a — > — oo as o? — > 1, and that a — > —1 as 
0? — > oo. Unfortunately the perturbation method (about d = 1) used by Collet et 
al gives only the behaviour of the largest eigenvalue 8. It will be observed from 
the graphs that all the eigenvalues appear to be smooth functions of o?, despite 
the fact that there are crossings in the exponents. Notice also that powers of 8 
appear in the scaling spectrum, and that these dominate the behaviour for d close 
to 1. (These powers of 8 can be seen as the rapidly increasing functions near o? = 
in Figure 3.) Conversely, the behaviour for large d is governed by a number of 
eigenvalues which approach each other in magnitude. In fact, these eigenvalues 
appear to occur in pairs of approximately equal magnitude but opposite sign. The 
numerical calculations suggest that the following limits exist, and take the very 
approximate numerical values given: 



lim 8il8 ^ 2 (23) 

lim 82/8 ^ -2 (24) 

lim ai/a ^ 2 (25) 

lim a2/a ^ —2 (26) 

4 Conclusion 

I have shown that all observable aspects of the period-doubling sequence in families 
of maps of the form A — | x | '^ are governed by the eigenvalues of the linearized period- 
doubling operator. This generalizes the classical Feigenbaum scaling. The same 
is no doubt true in the case of complex maps X — z"^ [m integer) in the case of 
n-tupling (n = 2, 3, 4, . . .). I have performed some preliminary calculations for the 
cases of period- doubling and period-tripling and m = 2, 3, . . . , 8 which confirm 
this. 
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